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The combination of Dirac physics and elasticity has been explored at length in graphene where 
the so-called ’’elastic gauge helds” have given rise to an entire new held of research and applications: 
Straintronics. The fact that these elastic helds couple to fermions as the electromagnetic held, implies 
that many electromagnetic responses will have elastic counterparts not explored before. In this work 
we will hrst show that the presence of elastic gauge helds will be the rule rather than the exception 
in most of the topologically non-trivial materials in two and three dimensions. In particular we 
will extract the elastic gauge helds associated to the recently observed Weyl semimetals, the ’’three 
dimensional graphene”. As it is known, quantum electrodynamics suhers from the chiral anomaly 
whose consequences have been recently explored in matter systems. We will show that, associated 
to the physics of the anomalies, and as a counterpart of the Hall conductivity, elastic materials 
will have a Hall viscosity in two and three dimensions with a coefficient orders of magnitude bigger 
than the previously studied response. The magnitude and generality of the new ehect will greatly 
improve the chances for the experimental observation of this topological, non dissipative response. 


I. INTRODUCTION 

In elasticity theory, Hooke’s law tells us that the stress 
tensor Tij applied on an elastic solid is proportional to 
the strain tensor Uij = {diUj + djUi)l2 as Tij = Xijwuir 
through a four-rank tensor Xijir- If the system is vis¬ 
coelastic the stress tensor is also proportional to the time 
derivative of Uij\ 

'^ij — Oijlr'^lr’) ( 1 ) 

through another four-rank tensor r]ijir. In general, this 
viscosity tensor possesses both symmetric and antisym¬ 
metric components under the permutation of pairs of in¬ 
dices. While the symmetric part is generally associated 
to dissipation and vanishes at zero temperature, the an¬ 
tisymmetric part arises when time reversal symmetry is 
brokerP. As usually happens with transport coefficients 
like the gyrotropic term of the permittivity in dielectrics, 
the antisymmetric part of the coefficients that are odd 
under time inversion are dissipation-less and they are not 
constrained to vanish at zero temperatures. In two di¬ 
mensions, isotropy, the symmetry under permutations of 
ij and Ir indices, and the aforementioned antisymmetry 
in permutations of pairs of indices allow for only one inde¬ 
pendent element of the antisymmetric part of the tensor 

Oijlr ' Vijlr ~ ^jlCr “k ^jrCl)-) whcrC Cij 

is the Levi-Civita tensor in two dimensions. 

The first context where this antisymmetric, 
dissipation-less coefficient was described was the 
quantum Hall effect (this is why this coefficient is called 
Hall viscosity)^. Being non-dissipative and appearing 
in quantum Hall effect systems, it is not strange that 
the Hall viscosity is associated to a certain Berry phase 
endowed with a topological meaning. In the case of 
the Hall conductivity, the off diagonal part of the 
conductivity tensor (Jxy is related to a non-zero Berry 


phase associated to the evolution of the phase of the 
wave function in the Brillouin zon^. In the case of 
the Hall viscosity, the parameter space is the space 
of the displacements of the atoms in the continuum 
limit, u{r). Since these displacements are considered 
adiabatic with respect to the electronic motion, one can 
consider a cyclic evolution of the wave function in this 
particular parameter space and define an associated 
Berry curvature of the wavefunction some energy 

level 

Also the Hall viscosity can be computed from a Kubo 
formulsP in terms of the stress-stress correlation func¬ 
tion Tjijir = f dte^^~^^(lTij(t),Tir(0)]} with the stress ten¬ 
sor defined as the variation of the Hamiltonian with the 
strain tensor T^j = After some mathematical ma¬ 

nipulations, a part of this Kubo formula appears to be 
proportional to the Berry curvature defined above. 

Besides, as it happens in all transport coefficients, the 
Hall viscosity appears in the effective action of some 
low energy degrees of freedom after integrating out the 
fermionic degrees of freedom. The effective energy func¬ 
tional for elastic displacements (phonons after quantiz¬ 
ing) for a system in two spatial dimensions consists on 
two terms, the standard elastic terms (possibly renor¬ 
malized by the effect of electrons) and the Hall viscosity 
term: 

a = ^ J(fr[XijirdiUjdlUr + rj^jirdiUjdlUr]- (2) 

Today there is a consensus about the meaning and the 
origin of the Hall viscosity. The question at this point is 
how to define the coupling between electrons and the elas¬ 
tic displacements Ui or better, to the strain tensor Uij and 
where to seek for it. The quantum Hall state is a topolog¬ 
ically ordered state so it is natural to extend the search to 
other topologically non-trivial quantum systems. Much 
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work has been done in this direction and Hall viscosities 
have been found in most of the topologically nontrivial 
condensed matter phases apart from the quantum Hall 
phas^ like the fractional quantum Hall phas^^, chiral 
superfluidJ® and superconductor^^, Chern insulators in 
presence of torsion in two dime nsionJ^^t^ and, more re¬ 
cently, in three dimension JISEH well. 

Precisely, the step of finding topologically non trivial 
systems from two to three spatial dimensions appears 
to be rather challenging. To have a non-zero Hall vis¬ 
cosity in three spatial dimensions breaking time rever¬ 
sal symmetry is not enough. Spatial isotropy or rota¬ 
tional invariance in the continuum limit must be bro¬ 
ken as welP. Up to now, only the torsional Hall viscos¬ 
ity defined in Weyl semimetalJI^ and three dimensional 
chiral superconductor^ has been discussed in the litera¬ 
ture. In these systems torsion is linked to the presence 
of dislocation^^ that define a preferred direction in the 
space. Also, in three dimensional topological insulators a 
Hall viscosity is allowed when the spatial isotropy is bro¬ 
ken by the presence of sample surfaces, and time reversal 
symmetry is broken by coupling magnetic elements to the 
two dimensional states appearing at these surfaceJ^. 

As we will discuss in detail later, the way how the 
elasticity couples to the low energy degrees of freedom 
(electrons or phase fluctuations in the case of superfiuids 
or superconductors) is crucial to determine if a topo¬ 
logically nontrivial phase displays Hall viscosity. Most 
of the studies available in the literature are based on 
a continuum formulation of the system and the elastic 
degrees of freedom. The philosophy is that, since elastic 
deformations can be viewed as geometrical deformations 
in the medium hosting the system, the excitations feel a 
distorted or curved space where to propagate. We can 
thus define an effective metric tensor related to elastic 
distortions and postulate that our system now develops 
its dynamics in a curved space. The general formalism of 
field theories in curved spaces and linear responses does 
the jotP^. Although powerful, this approach has some 
limitations in the scope of systems that can be treated. 
For instance, in the metric formalism (in absence of tor¬ 
sion), some non trivial topological phases with non-zero 
Hall conductivity have zero Hall viscosity. This is the 
case of Chern insulators in two dimensions and Weyl 
semimetals in three dimensions. On the other hand, 
only a small fraction of the literature devoted to the 
Hall viscosity in electronic systems treats the problem 
from the more conventional perspective in solid state 
physics through general electron-phonon coupling 
The goal in these approaches is to define an effective 
electron-metric coupling by taking the continuum limit 
from a lattice formulation, thus reaching similar success 
and suffering from the same limitations as in the former 
case. 

The present work offers a different alternative. We will 
show how the standard electron-phonon coupling in a 
class of electronic systems with an underlying lattice not 


only gives rise to a metric formulation in the continuum, 
but as we know from graphene, emergent vector fields 
arise in such effective description coupling to the elec¬ 
trons in a very similar way as the electromagnetic gauge 
field. These vector couplings between electrons and elas¬ 
tic deformations offer new results and enlarge the kind 
of systems where the Hall viscosity can in principle be 
measured. 


II. HALL VISCOSITY FROM HALL 
CONDUCTIVITY IN GRAPHENE 


Dirac materials is the generic name for electronic sys¬ 
tems whose low energy excitations are described by a 
Dirac Hamiltonian, graphene being the paradigmatic ex¬ 
ample. They exist in two and three spacial dimen¬ 
sions and have been the focus of attention in last years 
condensed matter. We will first show that anomalous 
Hall conductivity implies anomalous Hall viscosity in 
graphene and we will then extend the discussion to other 
Dirac materials. 

The low energy action around one of the Fermi points 
of graphene coupled to a U(l) electromagnetic back¬ 
ground field is 

^ ~ J — ivpi^di -h (3) 

A way to break time reversal invariance suggested by Hal¬ 
dane is to add complex next-to-nearest neighbours in the 
tight-binding on the honeycomb lattic^^^. In the contin¬ 
uum limit, these complex hoppings lead to a mass term 
with opposite sign at each Fermi point. After integrating 
out fermions it is easy to find the Chern Simons term as 
a part of the effective electromagnetic action: 

Scs = —sign{m) [(fxe^^^AidjAk, (4) 
TTC J 

after considering all the degeneracies (and in units of 
h = 1). The term @ is a consequence of the parity 
anomaly. Coupling the effective action @ to a back¬ 
ground electromagnetic source J^Ai allows to derive an 
anomalous quantum Hall effect: 

J* = — =- sign{m)e^^^djAk. (5) 

0 -^^2 ^ 


A very interesting aspect of graphene shared by many 
other 2D materials is the tight connection between lattice 
structure and electronic excitations. Within the simplest 
tight binding-elasticity approach, lattice deformations 
couple to the Dirac fermion as elastic U(l) vector fields 
t hat de pend linearly on the s train tensor uij = diUj+djUi 


(see section Methods V A[ ) 


Af = 


Af = 


^ {uii — U22) , 

- U12, 

a 


(6) 
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where p is the dimensionless Griineisen parameter esti¬ 
mated to be of order 2 in graphene, and a is the lattice 
constant. Since lattice deformations do not break time 
reversal symmetry, the elastic vector field couples with 
opposite signs to the two Dirac points and behaves as an 
axial vector field in a four dimensional representation. 

If we consider this elastic axial vector coupling to elec¬ 
trons instead of the electromagnetic field we have 


Sg = J d^x'ip — ivF'J^di -h -\- m'lp'ip, 


(7) 

It is obvious that the same line of arguments as before will 
give rise to a Chern-Simons action like 0 in terms of the 
elastic instead of the electromagnetic field. In particular 
we will get 


The same reasoning as done in the Haldane model can 
be applied to real graphene in a perpendicular magnetic 
field breaking time rever sal symmetry. The Hall conduc¬ 
tivity in graphene i^SHHIl = 4^2 l/ 2 )/ 7 r,n = 0 , 1 ,... 

By the same argument as before in the quantum limit 
where only the lowest Landau level is filled (n = 0), there 
is a contribution to the Hall viscosity coming from the 
elastic vector fields: 

VH ~ ( 12 ) 

The complete low energy action for deformed graphene 
in the elastic limit has been derived i]j 22 ] 

the presence 

of a magnetic field other terms coupling the elasticity to 
the Dirac fermions contribute to the Hall viscosity. In 
particular the following term 


= ‘^sign{m) J d^x (^AfAf - A^Af) , (8) 

substituting from eq. ^ we get in the Chern Simons 
action a term of the form 

Scs[y] = — 7 ^sign{m) / x[u22Ui2 - Uiiui2 
ira^ J 

+ U 12 U 1 I - Ui2U22]‘ (9) 


To this effective action we can now couple a source term 
for a time-dependent elastic deformation Su = UijT'^^ 
and compute the averaged value of the 11 component of 
the stress tensor 

(Jii) = = ^^sign(m)ui2. (10) 

ouli ira^ 

Hence we get the quantized Hall viscosity 

77 = — T^signym). ( 11 ) 

The coefficient is determined by the parity anomaly in 
( 2 - 1 - 1 ) dimensions, and the characteristic length associ¬ 
ated to the viscosity is given by the lattice spacing a. 
This new response is a genuine viscoelastic response: 
Elastic gauge fields inducing a Hall viscosity. 

The Hall viscosity described here is an intrinsic prop¬ 
erty of the Haldane model of graphene in the same sense 
as the quantum anomalous Hall effect (QAHE) of the 
original proposaP^. It has been arguecP^ that the intrin¬ 
sic Hall viscosity in the topological phase of the Haldane 
model is rjH = ^p/4, where p is the electronic density. 
When the Eermi level lies within the gap, no Hall viscos¬ 
ity appears in the Haldane model. 

Although the Haldane model has been realised experi¬ 
mentally in optical lattice^^^ it seems challenging to gen¬ 
erate it in graphene. The Haldane model is the model 
of the QAHE in graphene and the QAHE by magnetic 
proximity effect has been measured^, making graphene 
a suitable platform to potentially measure effects induced 
by the presence of the Hall viscosity. 


H = ^/3uij {'tp+aidj'ip + dj'ip^ai'ip ), (13) 

gives a contrib ution to the Hall viscosity that has been 
computed to b^^^ 


= 4/3^[|n|(|n| + l)/4 + l/8]eH/27r. (14) 


In the quantum limit (n = 0) we thus have 77 ^^ 
where Ib is the magnetic length. This result might seem 
a little bit puzzling at a first glance. We have stressed 
that in the absence of torsional couplings and magnetic 
fields, the Haldane model has no Hall viscosity unless one 
considers the elastic gauge fields. The reason is that the 
coupling in is a derivative coupling. When one con¬ 
sider both terms 0 and and find the effective theory 
for elasticity, the term coming from (13) has more than 
one space-time derivatives and it does not contribute to 
the Hall viscosity. When a magnetic field is taken into 
account, the partial derivative dj transforms into a co¬ 


variant derivative dj + ieAj (with B = \/ x A) and the 
magnetic length Ib = \JcjeB permeates through the cal¬ 
culations allowing for the result Taking graphene 

as an example, the ratio between the two contributions 
is 


Ph ^ 10 ^ 


(15) 


where a is the lattice constant and B is the magnitude 
of the magnetic field in Tesla. Eor a standard magnetic 
field of lOT, the contribution from the elastic gauge fields 
is three orders of magnitude larger than the one coming 
from (13). It means that if the Hall viscosity is measured 
in graphene under quantizing magnetic fields, it is the 
component discussed in the present work the one that 
will be measured. 


III. GENERALITY OF THE EFFECT 

The novel effect described in graphene is based on two 
main ingredients: The existence of elastic gauge fields. 
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FIG. 1: Two dimensional Dirac materials. The 

most representative example of the so-called Dirac 
materials in two spatial dimensions is graphene. In 
graphene, carbon atoms are arranged in a honeycomb 
lattice, (a). The highest occupied and the lowest empty 
bands touch at the corners of the Brillouin zone and the 
low energy band dispersion is well described by two 
species of massless Dirac fermions (c). M 0 S 2 is the 
representative of the family of transition metal 
dichalcogenides, (b). M 0 S 2 is a semiconductor with 
bands that admit a low energy effective description in 
terms of two species of massive Dirac fermions (d). 


and the non-trivial topology of the electronic bands of 
the system. This last ingredient is what allows to write 
down a Chern Simons effective action in 2D given in eq. 
0 which is the key issue of the approach (D will de¬ 
note the spacial dimensions of the (D+1) system). The 
non trivial topology can arise directly as a Berry phase 
associated to the Bloch bands of the crystal as in the 
Haldane modeP^, or from the topology of a general 2D 
electron gas in a pe rpen dicular magnetic field through 
the Hofstadter modeP^. In all cases the Hall conductiv¬ 
ity is associated to the topological Chern-Simons action 
depicted in eq. 0. As we have demonstrated, a Hall 
conductivity and elastic gauge fields inmediately implies 
the new Hall viscosity. We will now discuss the general¬ 
ity of the presence of elastic gauge fields in real systems 
and give some concrete examples. 

Although first derived in a tight binding formulation 
on the lattice, the generation of pseudogauge fields cou¬ 
pling shape deformations to electronic degrees of freedom 
in the low energy effective Hamiltonian of electron sys¬ 
tems, is a very general phenomenon. It also arises in a 
general construction of low energy effective actions based 
on a symmetry analysis as described irpfflUl^ xhe general 
conditions for a 2D lattice to support elastic gauge fields 
have been discussed iiP. The two essential ingredients 
are: The presence of (at least) two atoms per unit cell 
(this can arise from geometry as in graphene, from or¬ 
bital degrees of freedom, or by other mechanisms). Lack 


of inversion symmetry in the little group leaving a Fermi 
point invariant, ensures that the Fermi points sit at non¬ 
equivalent high symmetry points of the Brillouin zone. It 
is easy to see that the same structure giving rise to the 
Dirac Hamiltonian causes the minimal coupling to the 
vector fields associated to lattice deformations implying 
that the phenomenon can be generalized to other crystals 
in two and three dimensions. As it is clear from the anal¬ 
yses i] j^^ | 26 lthe symmetry group of the lattice dictates the 
precise form of the elastic vector fields. 

In 3D, an analysis of all crystallographic groups that 
sustain Weyl points in its low energy action in com¬ 
plete analogy to the case of graphene was first presented 
in ref.^. Since the Dirac structure in these crystals 
comes from the lattice, these having Weyl points in non¬ 
equivalent high symmetry points of the Brillouin zone 
will probably generate elastic vector field couplings un¬ 
der strain although a general analysis as the one done in 
2D is still lacking. 


A. 2D and Van der Waals systems 

New two dimensional materials of scientific and tech¬ 
nological interest are being synthesized in the “post¬ 
graphene” ersP^. They share many of the properties of 
graphene and often have large gaps making them more 
suitable for electronic applications. Single-layer and mul¬ 
tilayer transition metal dichalcogenides MX 2 {M = Mo, 
W and X = S, Se)P^ do support elastic gauge fields. 
Although having a noticeable band gap and a compli¬ 
cated orbital structure, M 0 S 2 and the related compounds 
based on the honeycomb lattice, support at low energy, 
the same elastic gauge fields as graphene given in eq. 
(|6|P^. These are also found in bilayer graphen^^ which 
has a quadratic dispersion and a Hall conductivity two 
times that of graphene. They will also occur in rombohe- 
dral multilayer graphene where the effective low energy 
model of the n layer compound has a Lifschitz dispersion 
relation uj ^ . Our mechanism for Hall viscosity will 

directly apply in these 2D (or effectively 2D) compounds. 

In these multilayer systems formed by a family of iso¬ 
lated metallic 2D lattice planes indexed by a primitive 
reciprocal lattice vector G^, the 3D Hall conductivity is 

_a6 ^ _2D^a6c/oO 

(Jh = —cfh e G ^. (16) 

Any 2D system having a quantum Hall effect (QHE) will 
exhibit, when stacked to form a 3D material, a 3DQHE 
if the inter-plane coupling is sufficiently weak, as it is 
the case of the so called Van der Waals (VdW) systems. 
When the 2D constituents support elastic gauge fields, 
a 3D Hall viscosity will arise whose coefficient is not re¬ 
lated to a quantum anomaly. The obvious example in 
this class is graphite whose 3D Hall conductivity has been 
described iiP^. In a homogeneous magnetic field perpen¬ 
dicular to the layers and when the Eermi level is in the 
3D gap, the Hall conductance is found to be (in units 
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FIG. 2: Band structure of representative Weyl 
semimetals. Typical band structure of a Weyl 
semimetal with two nodal points at different points in 
the momentum space. The vector connecting the nodal 
points is a key ingredient in the axial anomaly. 


of the quantum of conductance) (Jxy = where cq 

is the c-axis coupling constant of graphite estimated to 
be Co = 6 . 74 A. The corresponding 3D Hall viscosity will 
have an in-plane component 


^1112 — 2 


(2n + l)/3^ 
Tra^co 


(17) 


where a = 2.456Ais the in-plane lattice constant of 
graphene. Nowadays the family of VdW materials consti¬ 
tutes an intense area of research so we cannot limit our¬ 
selves to consider graphite or transition metal dichalco- 
genides. Hexagonal boron nitride (hBN) is another VdW 
system where these effects are potentially observable and 
even bulk black phosphorus whih has recently been found 
to have Weyl nodes under pressure}^ will also have the 
proposed elastic response. 


B. Anomalous Hall viscosity in a Weyl semimetal 

The derivation of the Hall viscosity from Hall con¬ 
ductivity can be extended to Weyl crystals in three 
spacial dimensionJ^. For electronic states in a three- 
dimensional periodic potential in a uniform magnetic 
field and when the Fermi level lies in an energy gap, the 
Hall conductivity adopts the fornP^ (Jij = ^^SijkGk^ 
where G is a vector in the reciprocal lattice of the peri- 

2 

odic potential. The Hall current is given by J = I^E x G. 
The presence of elastic gauge fields will give rise in these 
3D crystals, to an anomaly related Hall viscosity coeffi¬ 


cient rjH obtained from the 3D Chern Simons term 

Scs = J d^xGie^^'^^AfdkAf. (18) 

For a general 3D crystal breaking time reversal symme¬ 
try, the anomalous hall effect is characterized by a mo¬ 
mentum space vector v called the Chern vector. The 
anomalous hall conductivity is given by the same expres¬ 
sion: 


— 


27rc 


k 5 


(19) 


with the lattice vector replaced by z7. Our mechanism for 
the Hall viscosity will also apply in the three dimensional 
case by replacing the gauge fields in eq. (18) by the elastic 


gauge fields. The elastic Hall viscosity coefficient will be 


r]H 




( 20 ) 


where I is a characteristic length associated to the Chern 
vector. This is the basic mechanism at work. In (3+1) 
dimensions there are some subtleties concerning the coef¬ 
ficient of the Chern Simons action that will be discussed 
in detail in the Methods sec. IV B1 In what follows we 
will concrete these general expressions working out the 
anomalous Hall viscosity of Weyl semimetals. 

Weyl semimetals (WSM) are the “3D graphene”: 
They break time reversal symmetry and have pairs of 
Weyl nodes of opposite chirality connected by Fermi 
arcs. Their low energy description in terms of 3D Weyl 
fermions has given rise to an enormous amount of works 
related to quantum field theory physics: experimental 
realizati on of the axial anomaly, chiral magnetic effect, 
etc. (se d^^ l ^^ l and the references there in). To enhance 
the excitement, WSM physics has been r^orted recently 
in several compounds, particularly TaAP^. WSM whose 
Fermi surface crosses a given set i of Weyl nodes of chi¬ 
ralities located at positions Pi have an anomalous Hall 
conductivity given by eq. (19) with a Chern vectoJ^ 


z7 = Pi. Nevertheless the Hall viscosity of the in¬ 

trinsic material at zero temperature and chemical poten¬ 
tial is zercP^. By our proposed mechanism Weyl semimet¬ 
als will have an anomalous Hall viscosity at zero temper¬ 
ature and zero chemical potential given by the coefficient 
dictated precisely by the chiral anomaly supplemented 
with the elastic constant of the material. 

A detailed derivation of the elastic gauge fields for a 
lattice model of a Weyl semimetal is given in sec. Meth¬ 


ods V B The low energy effective action in the continuum 
limit around the Weyl points (±A) is given by {v = at) 

Hw{k) = (a{vkj_ ± A^l) + {v^ks ± Af)(T^) 

( 21 ) 

The elastic gauge fields are proportional to the strain ten¬ 
sor Uij with coefficients that depend on the parameters 
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of the model: 


Af = 
Af = 
Af 


-w? 

P\jb 3 -rv? 
r hi — r 

^^2 W 


'^31, 

'^32, 

■^33- 


( 22 ) 


Integrating out the fermions in the four components 
formulation where the field A couples as an axial field, 
we get through the electromagnetic chiral anomaly, the 
effective action 


■eff 




487r2 


63 — m ? 


J 


d^x {U 31 U 32 - U 32 U 31 ) (23) 


The numerical coefficient in (23) is explained in sec. 
Methods VB). The anomalous Hall viscosity (in the ab¬ 
sence of a magnetic field) is 


r]H = 


/3" 


1 

i3V ) ' 


247r^ a' 


(24) 


As it can be read from eq. (23) the new Hall viscos¬ 


ity obtained from the elastic gauge coupling is a coeffi¬ 
cient 773132 that points along the direction of the vector 
A breaking time reversal and rotation symmetries. Since 
in the perpendicular plane between the two nodal points 
the system is a Chern insulator (see Fig. [^, there will 
also be an in-plane component 771112 coming from the 
electron-phonon mechanism described in ref.^. 


IV. DISCUSSION AND WRAPPING UP 

In the present work, we have identified a new mech¬ 
anism to generate a Hall viscosity and a large number 
of new materials and systems that will display such ef¬ 
fect. As i t can be understood from the previous general 
discussiorPI^, it is natural to measure changes in the 
phonon structure due to the presence of 77 ^. In-plane 
phonon dispersion measurements can be performed by 
X-ray scatterin^^, Brillouin scatterin^^, and electron 
energy loss spectroscop}®l, all already done in the case 
of graphene and VdW materials. Also, since now we are 
also open to three dimensional electronic systems, three 
dimensional probes like neutron scattering can in prin¬ 
ciple be also used. The Hall viscosity can be related to 
changes in the phonon dispersions in presence of an ex¬ 
ternal magnetic field. In order to get some numbers we 
can take graphene as the paradigmatic example. Using 
standard valueJ^for the Lame coefficients A = 2 eVA“^ 
and /i = 10 eVA“^ a characteristic frequency associated 

2 

to the Hall viscosity can be defined, ujh = — 95 eV, 


which is a rather large frequency scale compared with the 
in-plane acoustic phonon frequencies in most systems (^ 
hundreds of meV). It intuitively means that the changes 
in the acoustic phonons will be hard to observe within 
current experimental resolutions at least in two dimen¬ 
sions. In spite of that, we believe that increasing the 
number of systems that display a Hall viscosity through 
elastic gauge fields, it is a matter of time to find an ex¬ 
perimental probe that could detect and measure 77 ^. 

We end up summarizing the main findings of this work 
and their physical implications: 

• The elastic gauge fields encoding the interaction 
of the electronic properties o with lattice deforma¬ 
tions of crystals are a very general phenomenon. 
First described in graphene, they have also been 
obtained in the low energy effective models of other 
2 D systems as M 0 S 2 but there were no examples in 
3D. We have deduced the presence and structure 
of the elastic gauge fields in a 3D model of Weyl 
semimetals (Sec. Methods). 

• In any compound supporting elastic gauge fields 
Hall conductivity implies Hall viscosity. This is 
particularly important in 3D where the precise form 
of the gravitational and chiral anomalies excludes 
the possibility of standard Hall viscosity in the ab¬ 
sence of torsion. For these systems with finite den¬ 
sity or temperature where a Hall viscosity exists, 
the contribution to the Hall viscosity from elastic 
gauge fields is several orders of magnitude bigger 
than that coming from phonons or metric deforma¬ 
tions. 

The lattice aspects of the viscoelastic response of topo¬ 
logical crystals have been recently analyzed iii^. Our 
new contribution to the Hall viscosity, although not ex¬ 
plicitly discussed, could certainly be worked out as a part 
of their general analysis. What we emphasize here is the 
analogy of the elastic gauge fields with the standard elec¬ 
tromagnetic partners which makes the connection with 
the anomaly related topological aspects very transpar¬ 
ent. 
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V. METHODS 


A. Elastic gauge fields in graphene 


The fact that the lattice deformations couple to the 
low energy electronic excitations of graphene in the form 
of elastic gauge fields has been explained in a number of 
articles and reviews. Nevertheless, being central to our 
work, we will repeat the derivation here taken mostly 
from ref.l^. Elastic deformations can be introduced in 
a tight binding calculation simply changing the hopping 
parameters. In the case of non equal hopping parameters 
7i = ^0 + ^2 5 oiie gets, near t he a K point, the effective 
low energy Hamiltonian 

H = vpip'^cFi {idi - Af) Ip, (25) 

where vp = (3/2)ta, t is the hoping parameter of the 
undeformed lattice and a is the lattice constant. The 
vector field A is 


Af 

Af 


-- (^2 + ^3 — 2ti) , 


(26) 


In the weakly deformed lattice, assuming that the atomic 
displacements u are small in comparison with the lattice 
constant a one has 

ti=t+ ^Si ■ {Ui - Uo), (27) 


where 6i are the nearest-neighbor vectors, uq is the dis¬ 
placement vector for the central atom, and 


Sint 
Sin a 


- 2 


(28) 


is the electron Griineisen parameter. The continuum 
limit (elasticity theory) with a displacement field u{r) 
is performed by making the substitution 


Ui-Uooc (ji • v) u (r^ , (29) 

from where we obtain the effective gauge fielcP^^^ 

Af = — {uii — U22) , 
a 

Af = -2 ^ui2. (30) 

a 


For the other valley, A', the sign of the vector potential 


(30) is opposite, in agreement with the requirement of 
time-reversal invariance. 

The gauge field (30) is proportional to the deformation 
tensor which is directly involved in the density of elastic 
energy. This means that, although the coupling to the 
fermions is gauge invariant, the problem as a whole is 
not. A pure gauge rotation by the gradient of a scalar 
function will change the kinetic term of the elastic part. 
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B. Emergent elastic vector fields in a Weyl 
semimetal. 

To illustrate how emergent vector fields associated to 
elasticity appear in a Weyl se mimeta l phase we can con¬ 
sider the following simple modeP^^of s-, and p-like elec¬ 
trons hopping in a cubic lattice and chirally coupled to an 
on-site constant vector field b. The parameters t, r, and 
m, represent, in a tight-binding description, the hopping 
matrix elements between s and p states, hopping between 
the same kind of states, and the difference of on-site en¬ 
ergies between s and p states, respectively. The vector 
field b breaks time reversal symmetry. Without loss of 
generality, we will choose the vector field b to point along 
the OZ direction. The tight binding Hamiltonian is: 

i^o = ^ {itaj - rP) Ci^j + (m -f 3r) ^ /3ci + h.c, 

i,j i 

(31) 

where i labels the position and j labels the six next 
nearest neighbors of length a in the cubic lattice. The 
matrices ai and (3 are the standard Dirac matrices. In 
the unstrained situation we will set all the hopping terms 
t equal for simplicity. 

We will focus on the parameter regime 0 > m > — 2r 
corresponding to a topological insulating phase. The long 
wavelength limit of this model around the T point {k = 0 ) 
is an isotropic massive Dirac system. Since Hq describes 
a massive Dirac model in the continuum limit in three 
spatial dimensions, in order to find a non vanishing Hall 
viscosity we need to break both time reversal symmetry 
and the effective rotational symmetry in the continuum. 
In the lattice model we can do both operations by includ¬ 
ing an on-site axial coupling between the fermion fields 
and a constant vector field b. Without loss of general¬ 
ity we will choose b to point along the third axis, with 
63 > m\ 


Hb = '^ bscf asjbCi. (32) 

i 


With the choice b^> m the spectrum of the total Hamil¬ 
tonian (31)+(32) consists in two bands crossing at two 
Fermi points in the Brillouin Zone and two bands at 
higher and lower energies, as it can be seen in Fig. <§■ 
To find the low energy effective model around these two 
Fermi points, we shall proceed in two steps: we will ex¬ 
pand the lattice model for small momenta compared with 
1/a, so after Fourier transforming (31 ) + ([3^ and intro¬ 
ducing the following approximations kja and 

cos{kja) 1 -one gets {v = at) 


Ho{k) = ^ {vttjkj +m^ + Ck = 

j 

= J2c+noik)ck. (33) 

j 


And now we project out the high energy bands. After 
a unitary transformation, the Hamiltonian matrix l-Lo{k) 


can be written as 

=(*'<'■ ‘•X + ih + vhn V 

^ ^ \ vk^cj^ uo-• + (63 - mjcTs ^ 

with k^ = (/ci,/c 2 ). The four-component wavefunction 
can be written in two-component blocks ( 0 ^, 7 /^). For 
energies E m ^ and momenta kj ^ (m + 63 )/^ 
we can write and the effective two-band 

model takes the form: 

Heff{k) =ip^ (va-k±-\ - -m? - v^kl)a3] tpk- 

\ m + D 3 / 

(35) 

From this effective Hamiltonian it is easy to see that 
the momenta where the two bands cross are X± = 

(0,0, Expanding now around these two 

points k X± E Sk, the final Hamiltonian takes the 
form of two massless three dimensional Dirac fermions 
'ip± separated 2A in the momentum space: 

Hw{bk) = {v<j ■ 5kji_ =F vsSksas) V’±.fe, (36) 

Now let us focus on the part of the original lattice 
Hamiltonian that depends on the strain. As discussed 
in the case of graphene, the strain tensor Uij enters in 
the tight binding approach through the change of the 
hopping parameters t when the lattice is distorted. The 
general recipe is (see also ref.l^) for ti, 

tai t{l - Puii)ai + tl3ui2a2 + t^uisas, (37) 

and the same for ti and ^ 2 , and 


- r{l - pUjj), 


(38) 


being p the corresponding Griineisen parameter of the 
model. Inserting these hopping changes in the original 
Hamiltonian (31), we can define the strained Hamiltonian 
as the sum of the original Hamiltonian Hq and the strain 
dependent part H[uij]. For our purposes we will first 
focus on these two parts of H[uij]\ 


5Hi[uij\ = tp'^cl[{uziai +Uz2a2 - ussas) sm{k3a)]ck, 
k 

(39) 

and 

5H2[uij] = r/3 ^ U 33 (l - cos(fc 3 a))c^^Cfc. (40) 

k 

Projecting out the high energy sector and expanding 
around the two nodal points A±, the strain dependent 
Hamiltonian part takes the form 

k i = 1.2 

- ^-^k^U33tp+,^a3'tp±,k- (41) 
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(a) (b) (c) 





FIG. 3: Evolution from a topologogical insulator to a Weyl semimetal. For 63 = 0 the spectrum consists in 
two pairs of degenerate bands due to time reversal symmetry, (a). When 0 < 63 < m the band degeneracy breaks 
down and a high energy sector differentiates from a low energy sector, but the system is still gapfull, (b). When 
63 > m, the low energy bands cross each other at two definite points in the Brillouin zone. At sufficiently low 
energies, the system consists in two pairs of Weyl fermions with opposite chirality, (c). 


We have found that, around the two nodal points, strain 
couples to the low energy electronic sector as a vector 
field: 


Af = Pyjbl- 
Af = P\Jhl- m?U32, 


A el 
^3 


^2 


m" 


-'^33- 


(42) 


We also note that this coupling is chiral, that is, the 
elastic vector field couples with opposite signs to the 
electronic excitations around the two Weyl nodes. What 
we have found is that, when the cubic symmetry in the 
original lattice model, or the continuous rotation symme¬ 
try in the effective model is broken by a constant vector 
6 , strain couples to electrons as a chiral vector field, sim¬ 
ilarly to what happens in graphene or other two dimen¬ 
sional systems. 

Integrating out fermions to get an effective theory for 
the elastic distortions will give rise, through the electro¬ 
magnetic chiral anomaly, to a term of the form 


rbid = - 


3 167r2 
1 

3 IGtt^ 


d^x (^AfAf - AfAf'^ = 


^7 

j d^x {U31U32 - U32U31 ), 


(43) 


and the Hall viscosity associated to the chiral anomaly 
mechanism is (in units h = 1) 


1 / 63 - \ I 


( 44 ) 


C. Chiral anomaly and Hall viscosity in Weyl 
semimetals 


Let us consider he action for a massless Dirac fermion 
coupled to a U(l) vector and to a constant axial vec¬ 
tor in (3-M) dimensions: 


7 


S= + eAf_^ + 


(45) 


At the classical level the action is invariant separately 
under vector and axial transformations: 7 /^ ^ exp[i(a]'0, 
7 /^ ^ exp [ 2 ^ 75 ] 7 /; what ensures the separate conservation 
of the vector and axial cur¬ 
rents: = 0 , = 0 . Also, since we have chosen 

to be constant, classically we can make a change in the 
chiral phase of the fermion fields to eliminate it from the 
action (45). As we know well, we cannot remove com¬ 


pletely the vector b^ from ou r the ory at the quantum 
level due to the chiral anomal^^EZl, that translates into 
the presence in the effective action for the electromag¬ 
netic field of the following term: 


ra[A] = ^J d^x{b,x^)e^^^^F^,F„0. 


(46) 


Integrating by parts and neglecting boundary terms, one 
can recognize in (46) the so-called Bardeen counter- 
ternP^. From the action (46), one easily obtains the 


anomalous current conservation laws at the quantum 


level: 9^ = 0 and 


1 nua0p p 
167r2^ 


ajS • 

Things become more involved when, besides the field 
A^{x) and the vector 5^, we consider adding to (45) an 
axial vector field that couples to the chiral current 
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(a) (b) 

FIG. 4: Anomalous Feynman diagrams. When both axial and vector fields are coupled to Weyl fermions, there 
are two diagrams that contribute to the axial anomaly. Imposing that the electromagnetic current must be 
conserved forces these two amplitudes to not be equal being the proportionality factor |. 


In this case the chiral current is not only coupled to and (46) is fixed by the requirement of gauge invariance 


a constant vector that leads to (46) after a chiral gauge 
transformation. In this case, the new effective action can 
be obtained by computi ng th e amplitudes diagrammati- 
cally depicted in Fig.(|4])f^^l^. Up to quadratic terms in 
and the effective action now has an extra term: 

ra[A^] = I d^x(b,xne^^^^UAx)fap{x), (47) 

where is the field strength associated to A^. The rel¬ 
ative factor of 1/3 in the coefficients of the terms in (47) 


= 0 . 

We are now ready to obtain the term in the effective 
action for elasticity corresponding to the Hall viscosity 
from (47) by identifying the vector field 6^ with (0, A) 
and A/as the elastic vector field in (42). According to 
our discussion on the elastic gauge fietes, since they are 
axial vector fields, and so is the vector 6, the coefficient of 
the effective action corresponds to the AAA vertex. This 
is the value that gives the final result for the anomalous 
Hall viscosity of the Weyl semimetal eq. (24). 













